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ABSTRACT 

We present rotation period measurements for 41 field M-dwarfs, all of which have masses inferred 
(from their parallaxes and 2MASS K-band magnitudes) to be between the hydrogen burning limit 
and 0.35 Mq, and thus should remain fully-convective throughout their lifetimes. We measure a wide 
range of rotation periods, from 0.28 days to 154 days, with the latter commensurate with the typical 
sensitivity limit of our observations. Using kinematics as a proxy for age, we find that the majority 
of objects likely to be thick disk or halo members (and hence, on average, older) rotate very slowly, 
with a median period of 92 days, compared to 0.7 days for those likely to be thin disk members 
(on average, younger), although there are still some rapid rotators in the thick disk sample. When 
combined with literature measurements for M-dwarfs, these results indicate an increase in spin-down 
times with decreasing stellar mass, in agreement with previous work, and that the spin-down time 
becomes comparable to the age of the thick disk sample below the fully-convective boundary. We 
additionally infer that the spin-down must remove a substantial amount of angular momentum once 
it begins in order to produce the slow rotators we observe in the thick disk candidates, suggesting 
that fully-convective M-dwarfs may still experience strong winds. 

Subject headings: stars: rotation - starspots - stars: low-mass, brown dwarfs - stars: evolution 

1. INTRODUCTION 

The rotational evolution of low-mass stars is predom- 
inantly governed by two competing processes. During 
the pre-main-sequence (PMS) phase, these stars are still 
collapsing, and thus, the moment of inertia decreases as 
a function of time. If angular momentum is conserved, 
the angular velocity must correspondingly increase. This 
spin-up, which persists until the star reaches the zero 
age main sequence (ZAMS), is counteracted by angular 
momentum losses, which are thought to be related to 
the star-disc i nteraction at early tim es (e.g. accretion- 
driven winds ; iMatt fc Pudritd 120051 or "disc locking"; 
Konigll 119911 : ICollier Cameron. Campbell fc Quaintrelll 
19951 k and stellar winds at late times, particularly after 
the star reaches the ZAMS. 

It is well-established that for solar-type stars, it is 
possible to reproduce the observed evolution reasonably 
well within this framework, taking as theoretical inputs 
a range of disc lifetimes in reasonable consistency with 
those observed in young clusters, a wind loss law involv- 
ing saturation of the angular mome ntum losses past some 
critical angular velocity aw , (e.g.. IStauffer fe Hartmannl 
Il987t iBarnes fe Sofia! fl996T ). and allowing the radiative 
core to decoupl e in angular velocity from the convective 
envelope (e.g.. iKrishnamurthi et aJJI1997t lAllainl IT9981 ; 
iDenissenkov et al1l2010D . A key feature of this solar- type 
evolution is that the spread in disc lifetimes gives rise to a 
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spread in rotation rates in young clusters, with the max- 
imal rotation rate being attained as the stars reach the 
ZAMS, after which they begin to spin down, and con- 
verge toward a narrow range of rotation rates. Around 
the age of the Hyades (625 ± 50 Myr from isochrone fit- 
ting; iPerryman et al.lll998| ) , the convergence is complete 
and all the stars follow a t 1 / 2 type spin-down thereafter, 
with rotation rate being a wel l-defined function o f mass 
and age for F, G and K stars (|Barnesll2003l l2007t) . 

Stellar winds depend on the magnetic dynamo that 
drives field generation, so the magnetic topologies of the 
stars are important in determining the angular momen- 
tum loss rates. For solar-type stars, this is thought to 
be an ail dynamo (where a and O refer to the "a ef- 
fect" , the twisting of the magnetic field lines caused 
by rotation, and the "tt effect", driven by differential 
rotation) operating at the interface between the radia- 
tive core and the convective envelope. For a fully- 
convective star, it is thought that this can no longcr 
operat e, and it has been suggested that a tu rbulent dy- 
namo ( Durnev. De Young fe Roxburgh! Il993) or a 2 dy- 
namo ( Radler et al.1 119901 ) may dominate. The former 
produces small-scale fields that would yield inefficient an- 
gular momentum losses through winds. 

In recent years, magnetic field measurements of stars 
on both sides of the fully-convective boundary have 
become available fe.g.. | Donati et al.l 1200 8; Mori n~et al.l 
120081 120101 : iReiners fc Basrill2009f ). These studies indi- 
cate there is indeed an abrupt change in field geometries 
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moving across the fully-convective boundary, with fully- 
convective objects storing more magnetic flux in large- 
scale fields than partially-convective objects. The effect 
of these changes on winds is not yet clear, but at any rate, 
we do not necessarily expect that the same wind formal- 
ism that works for solar type stars or partially-convective 
M-dwarfs should necessarily reproduce the observations 
for fully-convective M-dwarfs. 

While several authors have attempted to extend the 
analysis of rotational evolution to M-dwarfs, and par- 
ticularly to masse s below the fully-convect ive bound- 
ary (« 0.35 Mq; iChabrier fc Baraffd Il997r i. the lack 
of observations in this mass domain has made such 
an endeavor problematic. While observations in young 
open clusters have improved dramatically over the 
last decades, now yielding large samples of rota- 
tion periods below t he fully-convective boundary (e.g. 



1999; 
et al.l 



Herbst et all 120011: MakidoneTal] 



Stass un et al.l 
20041: iLamm et al.l 

2001 [2007; Sc holz. Eisloffel fc Munddl2009HCohen et al l 
2004 iLittlefair et al.ll2005l 120101; ICieza fc Baliberl 120061; 
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Irwin et al.ll2007bl . l2008af) . these extend typically only to 



a few 100 Myr, and at present only a handful of periods 
have been measured for older clusters, where the intrin- 
sic faintness of M-dwarfs makes determination of peri- 
ods difficult. Furthermore, such objects re ach the ZAMS 
at m uch later ages (a few hundred Myr; Bara ffe et al.l 
|1998|) than solar-type stars, so older clusters are needed 
to probe similar stages in the evolution. 

Additionally, the Sun has been used in all the studies of 
solar-type stars as a reference point to tie down the evo- 
lutionary models, particularly the normalization of the 
wind loss law, which relies on having data at late-times. 
For fully-convective stars, there is as yet no such object: 
a star with well known mass, age, and a robustly mea- 
sured rotation period (although Proxima comes closest to 
providing it; see later in this section). The best progress 
below the fully-convective boundary toward determin- 
ing the behavior at late times has been made through 
v sin ? observations of field st a rs (e.g., Delfoss e et al.l 
1998t iMohanty fc BasrH 120031; iReiners fc Basril 120081; 
Browning et al. 1 12010ft . While this information is invalu- 



able, and has led to important insights into the physical 
processes at play in these stars, it is important to note 
that v sin i observations are not sensitive to the slowest 
rotators, with typical spectral resolutions yielding a limit 
of 3 km s — 1 . Due to the small radii of M-dwarfs, this cor- 
responds to rather short rotation periods, e.g. 3.3 days 
for a 0.2 R star. 

Rotation period measurements of field M-dwarfs have 
provided a few clues as to what mi ght lie below the sensi - 
tivity limits of the vsini surveys. iBenedict et~all 0998) 
report a rotation period of 83 days for Proxima Centauri, 
a star which is highly likely to be fully-convective (they 
also report a weaker detection of a ~ 1 3 day periodicity 
in Barnard's star). iKiraga fc Stepienl (|2007l ) presented 
a comprehensive survey of rotation periods for field M- 
dwarfs, and whilst their sample is limited mostly to ob- 
jects above the fully-convective boundary (they present 
periods for two new objects below it), they were able to 
confirm a long period for Proxima, obtaining 82.5 days. 
Hart man et "all (|2009bfl also report a number of long rota- 
tion periods for field M-dwarfs from the HATNet transit 
survey. Finally, we see evidence for a comparably long 



rotation period in the tran siting exoplanet host GJ 1214 
(jCharbonneau et al.ll2009h . 

Th e MEarth transit s urvey (|Nutzman fc Charbonneaul 
120081; llrwin et al.ll2009D targets nearby (< 33 pc) north- 
ern hemisphere mid to late M-dwarfs to search for tran- 
siting super-Earth exoplanets in the habitable zones 
of their parent stars. The targets for this survey 
are selected to have inferred radii < 0.33 R© be- 
cause this is highly advantageo us for a transit survey 
(|Nutzman fc Charb onncau 2008) , and thus should all lie 
below the fully-convective boundary. Each star is typi- 
cally observed at moderate (20 minute) cadence on every 
clear night for a long period of time (~ 6 months with 
the current implementation of the survey). The data 
therefore possess favorable sampling to search for long 
rotation periods. 

2. OBSERVATIONS AND DATA REDUCTION 

MEarth is a targeted survey, and thus, has observa- 
tional properties rather different from the majority of 
transit surveys. We summarize the salient features in 
this section, particularly as they relate to the detection 
of continuous photometric modulations, rather than the 
discrete transit events for which the survey was designed. 

The data for the present work were gathered us- 
ing all 8 telescopes of the MEarth array, which is 
located within a single roll-off roof enclosure at the 
Fred Lawrence Whipple Observatory on Mount Hop- 
kins, Arizona. Data were from two full observing sea- 
sons, 2008/2009 and 2009/2010, which run from approx- 
imately mid-September to mid- July, where the observa- 
tory is shut down during the remaining « 2 month period 
corresponding to the summer monsoon season in south- 
ern Arizona. Each star was generally observed during 
only one of these two seasons, as our target list is cycled 
annually to increase the sample size searched for tran- 
sits. During each year of operations, hardware changes 
were minimal except for necessary repairs, and we used 
a fixed 715 nm long pass filter combined with a thinned, 
back illuminated 2048 x 2048 CCD on each telescope. 
The pixel scale is 0'.'76 per pixel, yielding a field-of-view 
of 26' on a side. 

Targets for the MEar t h sur vey were selected by 
iNutzman fc Charbonneaul ([2008) from a subsample of 
stars from the LSPM-North catalog (|Lepine fc Sharal 
l2005f ) with trigonometric parallaxes or spectroscopic or 
photo metric distanc e moduli indicating they are within 
33 pc (|Lepindl2005T) . 

For the purposes of this work, it is important to be able 
to, at a minimum, assign a reliable estimate of mass or 
spectral type to each target, and produce some estimate 
of the age. The majority of the MEarth targets are quite 
poorly characterized from existing data in the literature, 
often possessing only near-IR JHK magnitudes from the 
2MASS all-sky survey, photographic magnitudes from 
the Palomar Sky Surveys 1 , and proper motion informa- 
tion. We therefore elected to analyze only the « 1/3 of 
the sample with trigonometric parallaxes. There were 
273 such stars possessing more than one observation (in 
practice, more observations are of course needed to detect 



1 Note that although Lcpinc & Shara (2005) provide V-band 
magnitude estimates, the majority of these are from photographic 
plate measurements from the USNO-B1.0 catalog, and therefore 
have potentially large systematic uncertainties. 



Rotational evolution of fully-convective stars 



3 



a period; the smallest number of measurements for our 
selected rotation candidates was 207 for GJ 51, see also 
£|3.2[) . The combined parallax and K-band magnitude in- 
formation provides one of the best methods to estimat e 
the masses of sing le field M-dwarfs (jDelfosse et al.ll2000f ). 
and in combination with the proper motion information, 
can be used on a statistical basis to constrain the Galac- 
tic population to which the target belongs, and hence the 
age. Note that full kinematic information is not available 
for the majority of our targets as they lack radial veloc- 
ities, so this population assignment is necessarily quite 
crude at the present time. 

During the night, each telescope observes 20 — 30 fields, 
where the vast majority of fields contain only a single 
target M-dwarf. These are observed sequentially for the 
entire time they are above airmass 2, cycling around the 
list to yield a cadence per-field of ps 20 minutes. Expo- 
sure times are selected to yield a signal to noise ratio for 
the target star sufficient to detect a 2 R® planet tran- 
sit at 3cr per data p oint, for the as s umed stellar radius, 
which is based on the lDelfosse et all (|2000| ) K-band mass- 
absolute magnitude relation, a nd a po l ynom ial fit to the 
empirical mass-r adius data of iRibasI |2006) to convert 
mass to radius. 

Basic reductions and light curve generation were per- 
formed using an automa ted pipeline, based on that used 
by the Monitor project ()Irwin et al.ll2007al) . A number 
of instrument-specific refinements have been made, and 
these will be described in full in a forthcoming publica- 
tion (Berta et al., in preparation), along with our transit 
search procedures. There are, however, known to be two 
important systematic effects remaining after this stan- 
dard processing, which we will now describe as they are 
important for rotation period detection. 

MEarth uses German Equatorial Mounts, which ne- 
cessitate effectively rotating the telescope through 180° 
relative to the sky when crossing the meridian. Thus, 
each target samples two regions of the detector, one for 
negative hour angle and one for positive hour angle. Flat 
fielding errors are manifest as different base-line magni- 
tudes on each side of the meridian, so we solve for a 
"meridian offset" for each object to remove this effect. 

Secondly, we discovered correlations between the mea- 
sured differential magnitudes of the target M-dwarfs and 
weather parameters, specifically the relative humidity. 
This effect has been investigated in some detail and will 
be described by Burke et al. (in preparation). It results 
from a mismatch in spectral type between the target star 
and the comparison stars, where typically the compari- 
son stars are much bluer than the targets (usually, they 
are close to solar-type for the majority of our fields). 
Variations in the precipitable water vapor (PWV) con- 
tent of the atmosphere along the line of sight cause the 
strength of the telluric water vapor absorption in the 
MEarth bandpass to vary, which affects the photometry 
of the targets in a way which is not cor rected by stan- 
dard differential photometry procedures (|Angiond 119991 : 
iBailer- Jones fc Lamml 12003). This effect can reach sev- 
eral percent over the course of an observing season for 
the later M-dwarf spectral types. 

Since we lack sufficiently red comparison stars with 
good signal to noise ratios for the majority of our tar- 
get fields, we adopt an alternative method to derive the 
required photometric correction. We create a "common 



mode" light curve by taking the median of all of the M- 
dwarf differential magnitudes observed by all 8 MEarth 
telescopes in a half hour time bin (where the size was 
chosen to ensure all stars being observed at a given time 
are included). The averaging is necessary both to im- 
prove the signal to noise ratio, and to remove the effects 
of any real variability or transits. Since the precipitable 
water vapor usually varies slowly (the majority of the 
variations are from night to night) this is sufficient to re- 
move most of the effect. The mismatch in spectral type 
between the target and the comparison stars varies for 
each target, and therefore so does the amplitude of the 
effect, so in practice we subtract a scaled version of the 
"common mode", determining the scale factor for each 
object by standard least-squares methods. The scale fac- 
tor is generally found to be very well-correlated with the 
colors (and thus, presumably, spectral types) of the tar- 
gets. 

3. PERIOD DETECTION 
3.1. Method 

In order to properly account for the "common mode" 
systematic effect and the magnitude offsets between the 
t wo sides of th e merid ian, we adapt the method described 
in llrwin et al.l (|2006[ ). which uses least-squares fitting of 
sinusoids to the observed time-series m(t). 

We adopt as the null hypothesis a model of a constant 
(real) magnitude, modulated by the systematics correc- 
tions: 

TYi— + k c(t), h < 
m+ + k c(t),h > 



m (t) 



(1) 



where m_ and m + are separate (constant) baseline flux 
levels for the two sides of the meridian, and A: is a scale 
factor multiplying the "common mode" . h is hour an- 
gle, and c(t) is the "common mode" , determined from a 
time-binned median of all the M-dwarfs being observed 
by MEarth at any given moment. The scale factor k al- 
lows for the variable amplitude of the "common mode" 
component in each M-dwarf due to the differing spectral 
type mismatch of our targets and comparison stars. 

This null hypothesis is compared to the alternate hy- 
pothesis that the light curve contains a sinusoidal mod- 
ulation, again modulated by the systematics corrections, 
of the form: 

, s _ J ui— + k c(t) + asm(u)t + <j)) , h < /0 -, 
mi(t > ~ \m+ + k c(t) + a sin(wi + (j)),h>0 W 

where a and <fi represent the semi-amplitude and phase of 
the sinusoid, and w = 2tt / P is the angular frequency cor- 
responding to rotation period P. In order to determine 
uj, we fit this model using standard linear least-squares, 
at discrete values of uj sampled on a uniform grid in fre- 
quency from 0.1 to 1000 days. By comparing the best- 
fitting x 2 values for the two hypotheses as a function of 
frequency, we derive a "least-squares periodogram" that 
accounts for the effect of the systematics. 

Due to the small sam ple size, we omitte d the cut in 
X 2 improvement used by llrwin et al.1 (|2006f l. and simply 
subjected all of the light curves to visual inspection, to 
define the final sample of periodic variables, although in 
practice all of the objects s e lected pass the Axl > 0.4 cri- 
terion used by llrwin eTaTI (l2006h . 41 light curves passed 
this selection, where the remainder were consistent with 
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no detectable variation, had excessive systematics, or had 
insufficient data to determine a period. A number of ob- 
jects appeared to have monotonic trends in time, but it is 
not clear for many of these if they are due to systematics 
or to variability at present. 

3.2. Sample properties 

Figure Q] summarizes the overall properties of the sam- 
ple of 273 stars on which the period search and selection 
of rotation candidates was performed. As discussed in 
$2] all of these stars have trigonometric parallaxes, and 
Figured] shows that all of the objects selected as rotation 
candidates also have more than 200 data points taken on 
> 10 nights (for all but one object on > 28 nights), span- 
ning at least 120 days. 

It is important to note that the selection we have per- 
formed does not remove close binaries from the sample. 
Tidal effects will modify the rotation rates of the com- 
ponents of such binaries, by transferring angular mo- 
mentum between the stellar spin and the binary orbit. 
The dependence of tides on binary semimajor axis is 
very strong (e.g. IZahdfl977t ) so the objects showing the 
largest tidal effects will be in the shortest period systems. 
It is therefore possible that some of the objects in the 
present sample have had their rotation rates altered by 
a binary companion, with the objects rotating at shorter 
periods being more likely to have been affected. 

We also note that the use of K-band absolute magni- 
tude to estimate masses means that these masses will be 
overestimated for any unresolved binary or multiple star 
systems due to the extra light from the companion(s). 

These issues are best resolved by performing spectro- 
scopic follow-up at high resolution, to search for double- 
lined objects and for radial velocity variability in order 
to identify any spectroscopic binaries. This has not yet 
been done for the present sample, and it is important to 
bear the caveats discussed in this section in mind in the 
interpretation of the rotation periods we measure. 

3.3. Effect of the "common mode" correction 

Since we are solving for the "common mode" ampli- 
tude simultaneously with the sinusoidal variability, it is 
important to evaluate the effect of this correction on our 
period sensitivity. 

Figures [5] and [3] show periodograms of the "common 
mode" light curve, calculated in the same way as those 
we use for period detection by least-squares fitting of 
sinusoids. The dominant power in the "common mode" 
is on long timescales, and at the corresponding 1 day -1 
alias frequencies, with the highest peaks for the two years 
of observations corresponding to periods of 25.1 days for 
2008/2009 and 14.5 days for 2009/2010. 

The presence of power at such frequencies could af- 
fect our sensitivity to rotation periods in this range. We 
therefore proceed in the next section to simulate the full 
end-to-end period detection process, including the ef- 
fect of the "common mode" and of the visual inspection 
stages, to evaluate the survey sensitivity. 

3.4. Simulations 

In order to evaluate the sensitivity in period and ampli- 
tude, simul ations were perform ed following the method 
detailed in llrwin et al.l (|2006f) , injecting sinusoids with 
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Fig. 1. — Summary of sample properties. In each panel, the open 
histogram shows the full sample of 273 stars with more than one 
observation, and the solid histogram shows the 41 stars for which 
we detect periods. The panels show the number of observations 
for each target (top), number of nights during which data were 
gathered (center), and the range of Julian dates spanned by the 
observations (bottom). 

periods from 0.1 to 200 days following a uniform distri- 
bution in log period into only the light curves rejected 
in the previous visual inspection stages to reduce con- 
tamination by any real variability. Two semi-amplitudes 
were simulated, 0.005 and 0.01 mag, corresponding to the 
range of typical amplitudes of our rotation candidates. 

The results are summarized in Figures |4] and [5] for the 
0.005 and 0.01 mag semi-amplitudes, respectively. We 
also show plots of the recovered period versus injected 
period, in Figures |5] and [71 respectively. Whilst the over- 
all completeness of rotation period detections is quite low 
(around 50 — 60% in most period bins), there is no clear 
bias except in the longest-period bin, where the com- 
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Fig. 2. — Periodogram of the "common mode" light curve for 
the first year of observations. Plotted as the square root of the 
reduced x 2 of the sinusoidal fit as a function of frequency, with 
the vertical axis inverted such that higher peaks correspond to 
more significant detections as for a conventional periodogram. The 
period of the highest peak is 25.1 days. When computing the x 2 
values we used the scatter in each time-bin to compute the weights; 
this appears to significantly underestimate the scatter between the 
bins, hence the large values of the reduced x 2 - The peaks around 
1 day - appear to be caused by aliasing, and probably do not 
represent real variations in the PWV on these timescales. This has 
been confirmed by analyzing contemporaneous GPS-based PWV 
measurements (which are not restricted to hours of darkness) from 
Flagstaff, Arizona. 
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Fig. 3. — As Figure [2] only for the second year of observations. 
The period of the highest peak is 14.5 days. 

pleteness drops to 25 — 30%. This is expected due to the 
limited survey duration. 

The reliability and contamination histograms, and Fig- 
ures [H] and [7| provide an indication of the reliability of 
period recovery in the cases where a significant modula- 
tion was detected. For 0.01 mag, these statistics indicate 
very good period recovery, with reliability (measuring 
the fraction of objects detected with the correct period) 
above 80% in all period bins, and contamination below 
20%. For 0.005 mag the performance is substantially 
worse, with a large drop in reliability around 1 day and 
in the longest-period bin. The contamination statistic 
shows a similar effect. By examining Figure H2 it is clear 
that some of the scattering of objects in and out of the 
1 day bin is due to aliasing, and this diagram also indi- 
cates that many of the "incorrect periods" contributing 
to the poor performance for the longest periods simply 
have larger period errors than our 10% threshold. This 
is probably the reason for the apparently small change 
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Fig. 4. — Results of the simulations for 0.005 mag semi-amplitude 
expressed as percentages, plotted as a function of period. Top 
panel: completeness as a function of real (input) period. Center 
panel: Reliability of period determination, plotted as the fraction 
of objects with a given true period, detected with the correct period 
(defined as differing by < 10% from the true period). Bottom 
panel: Contamination, plotted as the fraction of objects with a 
given detected period, having a true period differing by > 10% 
from the detected value. 

in the completeness statistic between the two bins: this 
merely counts detections, without regard to whether the 
period was correctly determined. 

While it is important to perform the simulation on M- 
dwarf target stars to account for the systematics, it is 
expected that in reality, many or even the majority of 
these should show real, astrophysical variability at some 
level. It is therefore likely that we have injected our sim- 
ulated signals into objects which also have astrophysical 
signals. For the 0.01 mag amplitude, this does not ap- 
pear to be a serious problem, with the injected signal 
generally overwhelming anything already present. How- 
ever, this could contribute to explaining the apparently 
high contamination and the lack of a significant drop in 
completeness in the 0.005 mag amplitude sample. 

In order to investigate these effects, and in particular 
to also test the influence of systematics on the detec- 
tion rate, we have performed an additional set of sim- 
ulations for 0.005 mag semiamplitude but using white, 
Gaussian noise of standard deviation set by the estimated 
observational errors, rather than using the observed light 
curves. This procedure eliminates systematics ("corre- 
lated noise") and any existing astrophysical variability 
but should otherwise have a comparable noise level to 
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Fig. 5. — Panels as Figure[4]except for 0.01 mag semi-amplitude. 




Input period (days) 

Fig. 6. — Recovered period as a function of actual period for 
our simulations at 0.005 mag semi-amplitude. Crosses represent 
objects where the correct period was recovered to within 10%, and 
open circles show objects failing this criterion. The diagonal line 
indicates equal periods, and the curves show the loci of the 1 day -1 
aliases. 




Input period (days) 

Fig. 7. — Recovered period as a function of actual period for 
our simulations at 0.01 mag semi-amplitude. Symbols and lines as 
Figure [6] The arrow at the top of the diagram indicates a single 
point which fell outside the plotted range. 

the real data. The results are shown in Figure [5J This 
indicates that the low completeness seen in Figures 2] and 
[5] at short periods results from a roughly equal contribu- 
tion of systematics (or contamination from variability) 
and other effects (sampling and noise). The completeness 
at long periods is slightly improved, but as expected, is 
still lower than at short periods, remaining at the ~ 50% 
level. This is most likely due to the sampling. The con- 
tamination and reliability histograms show a substantial 
improvement, with close to ideal period recovery, which 
suggests these effects are almost entirely due to system- 
atics (or contamination from variability). 

4. RESULTS 

Our results are presented in Table [TJ and the phase- 
folded light curves of all 41 objects are shown in Figure 

El 

It is evident from Figure [S] that the light curves for 
several objects (e.g. LHS 1068, G 32-37A, GJ 170, LHS 
543a) are not perfectly sinusoidal. This is indeed ex- 
pected if the modulations arise from stellar spots. It is 
also important to note that there exist pathological spot 
distributions (e.g. two identical active longitudes spaced 
by 180°) which can lead to doubling (or in the general 
case, multiplication by an integer) of the frequency, caus- 
ing us to mis-estimate true the rotation period of the star. 
One way to check for this is to compare the rotation pe- 
riods with vsmi measurements. 

We have searched for literature measurements of ra- 
dial and rotational velocities for our targets, finding 
ysini measurements for 8 objects (iDelfosse et all 119981 : 
iMohantv fc Basril2003UBrowning et al.ll2010[ ). These are 
summarized in Table [2] which includes radial velocities 
and kine matic population as signments for two additional 
objects (jMontes et al.ll2001h . 

Figure [TU] shows a comparison of these v sin i measure- 
ments with the rotation velocities inferred from our pe- 
riod measurements and radius estimates. These are gen- 
erally in good agreement, with the exception of one ob- 
ject, LP412-31. This is one of the latest-type objects in 
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TABLE 2 

Rotational velocities, radial velocities and kinematic populations 
from the literature for our targets 



CD O 



Name 


v sin i 
(km s —1 ) 


a[y sin?) 
(km s _1 ) 


^rad 

(km s _1 ) 


Population 


h 

Reference 


GJ51 






4.4 


YD 


1 


GJ1057 


< 2.2 




27 


OD 


2 


LP412-31 


12.0 


2.0 


48.8 


YD 


3 


G099-049 


7.39 


0.8 


30 


YD 


2 


LHS1885 


< 3.7 




16 


YO 


2 


GJ285 


4.6 


0.4 


26.5 


YD 


4 


GJ1156 


9.2 


1.9 


4 


YD 


2 


GJ493.1 


16.8 


2.1 


-27 


YD 


2 


GJ791.2 


32.0 


2.0 


-29 


YD 


2 


GJ856A 






-24 


YD? 


1 



Sec Leggett ( 1992 



b (11 iMontes et al.1 (p0 0l : (21 iDelfosse et al.l flUl); (3) 
iMohantv fc B asri (|200l : (4) iBrowning et~all pOlOh 





Period (days) 

Fig. 8. — Panels as Figure [4] except simulating light curves with 
white, Gaussian noise rather than using real data. 



our sample (|Mohantv &: Basrill2003l state a spectral type 
of M8.0), so we suspect the discrepancy may be caused 
by an error in the assumed radius when calculating v m t, 
especially no ting the lack of objects at these late spectral 
types in the [Ribas (2006) sample we used to derive the 
empirical mass-radius relation. 

4.1. Kinematic population and age assignment 

It is important for the analysis of rotational evolution 
(performed later in this section) to be able to assign ap- 



proximate ages to the MEarth sample. This is a notori- 
ously difficult problem for field stars. For the purposes 
of the present work, we use the available kinematic in- 
formation (5 of the 6 phase-space dimensions, where we 
lack radial velocities for the majority of the targets) to 
statistically assign our targets to the Galactic thin disk 
or thick disk / halo populations. 

Previous studies ha ve used tangential velocity (e.g. 
iReiners fc Basril [20081 ). but this does not take into ac- 
count the bulk motions of the various Galactic popula- 
tions along each line of sig ht through the Ga laxy. In- 
stead, we use the model of iDhital et al.1 (|2010[ ) to com- 
pute the predicted mean proper motions and their dis- 
persion for the Galactic thin disk at the position of each 
M-dwarf, and compare to the measured proper motions. 
Stars within < la of the thin disk prediction were la- 
beled "thin" (likely thin disk members), stars at > 3a 
were labeled "thick" (corresponding to the thick disk / 
halo) and are likely to be old, and the intermediate cases 
were labeled "mid" . These population assignments are 
reported in Table [TJ 

We have assigned a 0.5 — 3 Gyr age to the "thin" 
stars, and a 7 — 13 Gyr age to the "thick" stars (e.g., 
iFeltzing fc Bensbvl 120081) . The "mid" sample in reality 
contains a mixture of thin and thick disk populations, so 
we assign it an intermediate age of 3 — 7 Gyr. 

4.2. Rotation period distribution 

Figure [TT] shows rotation period plotted as a function 
of stellar mass, where the symbols indicate our kinematic 
population assignments from ij4.ll This diagram shows 
no clear trend of rotation period with stellar mass. The 
population appears to be in two main clumps in the dia- 
gram, a population of very rapidly rotating objects with 
periods of ~ 0.2—10 days, and a population of very slowly 
rotating objects with periods of 30 days to the sensitivity 
limit of the observations at approximately 150 days. 

The "thin" sample falls entirely within the rapidly ro- 
tating clump of objects, whereas the majority of the 
"thick" stars fall in the slowly-rotating clump. The sit- 
uation is intermediate between these extremes for the 
"mid" sample. Given the likely ages of these populations, 
this indicates that the older objects in the sample are on 
average rotating more slowly. Such a trend is a natu- 
ral consequence of the expectation that these stars spin 
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GJ1006B P=104d 
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Fig. 9. — Phase-folded light curves for our 41 rotation candidates. Differential maj 
In these diagrams, we plot <j>/2Tr on the horizontal axis, where <j> is defined in Eq. 
"common mode correction" described in the text. 



mitudes are plotted, in the MEarth instrumental system. 
J2J, and have subtracted the "meridian correction" and 
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TABLE 3 

List of references for rotation periods from the 
literature. 



V rot ( km / s ) 

Fig. 10. — vsini from spectroscopic measurements plotted as 
a function of rotation velocity from our period measurements, as- 
suming the stellar radii given in Table [T] Points with arrows show 
the two vsini limits. Open symbols denote young/old (YO) and 
old disc (OD) objects, and filled symbols denote young disc (YD). 
The gray shaded region indicates the part of the diagram where 
vsini > v TO t, and should not be occupied. Note that the stellar 
radii (and thus the estimated rotation velocities) also have signif- 
icant uncertainties, but these are difficult to quantify, so we have 
not plotted horizontal error bars on the diagram. These uncer- 
tainties are likely to amount to at least 10%, and probably greater 
give n the uncertainties in the effects of activity on M-dwarf radii 
(e.g. IChabrier efaHl200?D . 



Mass (M Q ) 

Fig. 11. — Rotation period as a function of stellar mass for our 
sample. The symbols indicate our kinematic population assign- 
ments: open circles indicate "thin" objects (likely thin disk and 
hence younger), gray filled circles "mid" objects, and black filled 
circles indicate "thick" objects (likely thick disk or halo and hence 
older) . 

down due to winds as they age, and the lack of slowly ro- 
tating stars in the "thin" sample compared to the abun- 
dance of them in the "thick" sample further indicates 
that this spin-down takes place during ages intermedi- 
ate between these two samples for the majority of the 
M-dwarfs. These c onclusions confirm the findings of the 
studies (e.g. iDelfosse et aT1ll998t E einers fc Basril 



Cluster 


Age 
(Myr) 
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1 


Hcrbst et al. (2001. 2002) 
Stassun et al. (1999) 


NGC 2264 


2 


Lamm et al. (2005) 
Makidon ct al. (2004) 


NGC 2362 


5 


Irwin ct al. (2008b) 


IC 2391 
IC 2602 
NGC 2547 


30 
40 


Patten & Simon ( 1996) 
Barnes et al (1999T~ 
Irwin ct al. (2008a) 


Pleiades 


100 


Hartman ct al. (2009a) 
Terndrup et al. (1599) 
Scholz k Eisloffcl (2004b) 


NGC 2516 
M35 
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Irwin ct al. (2007b) 
Mcibom ct al. (2009) 
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Praesepe 
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Radick ct al. ( 1987) 
Prosser et al. ( 1995) 
Scholz & Eisloffcl (2007) 
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Baliunas, Sokoloff & Soon (1996) 
Barnes (2001. 2003) 
Kiraea Oteoich (2007) 
Benedict et al. ( 1998) 









I2001 . 

Using a two-sided Kolmogorov-Smirnov test, we have 
verified that the difference between the "thin" and 
"thick" periods is statistically significant, with a prob- 
ability of 4 x 10 ~ 3 that these are drawn from the same 
parent population. The corresponding probabilities for 
the "thin" and "mid" , and the "mid" and "thick" pairs, 
are 0.13 and 0.06, respectively, neither of which we con- 
sider to be statistically significant. This is not surprising 
as the sample sizes are quite small and the "mid" popu- 
lation likely consists of a mixture of the other two pop- 
ulations. The median periods are 0.7, 7.8, and 92 days 
for the "thin" , "mid" and "thick" samples, respectively. 

5. DISCUSSION 
5.1. Morphological comparison with open cluster data 

In Figure rTJ] (see Table [3] for a list of data sources) we 
show Figure [TT] in the wider context of rotation period 
distributions from the literature for younger stars in open 
clusters, and for more massive stars in the field. 

For the MEarth sample, we have assumed the "thin" 
population and "young disc" are equivalent, and likewise 
for "thick" and "old disc" . The "mid" objects are plotted 
in both panels using open symbols as it is not clear to 
which population these should be assigned, and it is likely 
the "mid" sample contains members of both, with the 
slower rotating objects being more likely to belong to 
the "old disc" population and vice versa. 

Figure [12] shows that the pre-main-sequence behavior 
of M-dwarfs is markedly different from that of solar-type 
stars. The open cluster data indicate that the M-dwarfs 
all appear to spin up rapidly, and reach very rapid ro- 
tation rates at 100 — 200 Myr age, which corresponds 
to the time at which the most massive objects in the 
< 0.35 M domain reach the ZAMS. Furthermore, the 
slowest rotation period seen at any given mass is itself 
a very strong function of mass, declining moving to less 
massive objects. The similar morphologies in the dia- 
grams, especially from NGC 2264 to NGC 2516 (after 
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Fig. 12. — Compilation of rotation periods for low-mass stars from the literature, including the results from the present study. Plotted 
in each panel is rotation period as a function of stellar mass for a single cluster, or several clusters that are close to coeval. In the case 
of the field star samples, we divided them into two bins for the "young disc" and "old disc" stars, respectively, and color all the literature 
data gray to distinguish them from the results of the present study (in black). We plot the "thin" population in the "young" panel and 
the "thick" population in the "old" panel, both as solid symbol s, with the "mid" p opulation shown in both panels using open symbols (sec 
text) . Measurements for field sta rs are from the compi lation by Barnes (2001, 2003) of the rotation periods for the Mount Wilson stars from 
Baliunas. Sokoloff & Soon (1996), down to pb 0.7 Mg, Kiraga & Stepicri (2007) for early-M dwarfs, and the present study for later objects. 
The appropriate references for each panel are given in Table [3] A ll the masses for the cluster samples were computed using the I-band 
luminosities of the sources and the models of [Baraffc ct a f] (119981), assuming values of the age, distance modulus and reddening for the 
clusters taken from the literature. For the field stars, the lDelfosse et al. (2000) mass-absolute magnitude relations w ere used for M- dwarfs, 
and for all higher-mass stars, we converted t he observed B — V c olors to effective temperature using Table A5 from Kcnyon & Hartmann 
(1995), and then to mass using the models of [Baraffc ct al. (19§3). 



accounting for the expected spin-up due to contraction, 
which corresponds simply to a vertical translation in this 
logarithmic plot), indicate there is relatively little rota- 
tion rate dependence in the PMS evolution. 

Betwe en the Pleiades age 100 Myr from isochronc 
fitting; iMevnet, Mermilliod fc Maederl 119931) clusters 
and the much older Praesepe cluster, the available data 
indicate there is essentially no evolution of the rotation 
period, so it appears that little angular momentum is 
lost over this phase of the evolution, corresponding to 
the arrival at the ZAMS and the early main sequence for 
these stars. It is, however, important to note that the 
available samples at the Hyades/Praesepe age are very 
limited at the present time, so we must caution drawing 
strong conclusions from them. We shall return to this 
point later in the discussion. 

While the M-dwarfs do not appear to have spun down 
at the age of the Hyades/Praesepe, this process is clearly 
well underway for the earlier M-dwarfs by the "young 



disc" age, and apparently complete by the "old disc" age 
above approximately 0.4 M Q . Below this mass, there is 
a dramatic change moving from the "young disc" sample 
to the "old disc" sample (especially recalling the uncer- 
tainty as to which bin the "mid" objects, plotted as open 
symbols, should be assigned), where many slowly rotat- 
ing stars are now found, but there is still a wide range of 
rotation periods here, in contrast to the higher masses, 
spanning about 2 orders of magnitude. 

It is interesting, and suggestive, to compare the rota- 
tion period distribution in the lower right pan el of Figure 
12l w ith the activity fraction (e.g. Fig. 3 of We st et all 
2008) . The fraction of active field M-dwarfs shows a steep 
rise around M5 spectral types. It is at approximately this 
point moving down in mass where the rapid rotation phe- 
nomenon appears at the oldest ages. Given the age range 
displayed in the diagram, this seems most likely to be an 
evolutionary effect: the time at which the stars spin down 
increases with decreasing mass, becoming comparable to 
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Fig , 13. — Evolution o f the moment of inertia from the mod- 
els of|Baraffe ct al. (1998) for the radiative core (lower lines) and 
convective envelope (upper lines) for masses of 1.0 (solid line), 0.8 
(dashed line), 0.6 (single dot dashed line), 0.4 (dotted line), and 
0.2 Mq (triple dot dash line). Note that the 0.2 Mq star does not 
form a radiative core so only one line is shown. 

the sample age below the fully-convective boundary. 

The existence of a wide range in rotation rates, and 
the appearance of a substantial number of slow rotators 
below the fully-convective boundary moving between our 
"thin" and "thick" bins (although some are also evident 
in the "mid" bin) further suggests that any spin-down 
occurs quite rapidly in the interval between the nominal 
ages of these samples for the majority of fully-convective 
objects, i.e. around 5 Gyr in age. 

These timescale s are in good agreement with those of 
IWest et all (|2008| ). e.g. their Fig. 10, where they find 
activity lifetimes of ~ 7 Gyr below the fully-convective 
boundary, and < 2 Gyr for the partially-radiative objects 
above it. We shall explore possibilities for generating 
such a spin-down in terms of wind models in £15.21 

We note some morphological features in the diagram 
for the "old disc" sample. The hi g hest m ass stars show 
the familiar break in the Kraft (1967) curve around 
1.2 M©, which is thought to be due to the disappear- 
ance of the surface convective zones moving to higher 
masses. Stars without surface convective zones appear 
to experience no (or very little) wind losses on the main 
sequence and thus remain at very rapid rotation rates 
until they start to evolve off the main sequence. 

There is another discontinuity in the diagrams around 
0.5 - 0.6 M , which is highly evident in the NGC 2516 
distribution, and has been notic ed and discussed by many 
authors. We note here that the lKiraga &: Stepieril ()2007f ) 
sample appears to also show a discontinuity at this mass 
for field M-dwarfs, where the rotation period has rel- 
atively little mass dependence immediately above this 
mass, and appears to increase with decreasing mass be- 
low it. 

5.2. Rotational evolution models 

Figure [14] gives an alternative presentation of the avail- 
able rotation period data for objects below 0.35 M©, in a 
way where the evolution is more explicit, at the expense 
of the mass information seen in the period versus mass 
plots. It is important to bear in mind the strong mass de- 
pendence, particularly in the open cluster samples, when 
trying to interpret this figure. We have assigned ages to 
the MEarth sample based on the discussion in $4.11 not- 



o 

~ : T di3CS | DW = 4Myr w sQt = 3.'l 3w - 

' T di 3C f Q5 t = 16M y r ' K wsto w = 1-1e+47 | 

■ ■ + + K w'fa 5 t = 1.2e+45 J ~ 




Proximo 

N2516 Praesepe 

ONC 2264 2362 N2547 Pleiades Thin Mid Thick ' 

Barnard 



I 10 100 1000 10 4 

Age (Myr) 

Fig. 14. — Rotational angular velocity u) plotted as a function 
of time for stars with masses 0.1 < AZ/Mq < 0.35. Crosses show 
the rotation period data such that each cluster (or field sample) 
collapses into a vertical stripe on the diagram, and short horizontal 
lines show the 10th and 90th percentiles of lj, used to characterize 
the slow and fast rotators respectively. Plotted are the ONC, NGC 
2264, NGC 2362, NGC 2547, the Pleiades, NGC 2516, Praesepe, 
and the MEarth data. We also show t he po sitions of Proxima 
and Barnard's star from [Benedict ct al. (1998) as the large, bold 
crosses. We have assigned ages to the field stars for the purposes 
of plotting them on this diagram, assuming 6 Gyr for Proxima 
(the center of the range of ages published for the a Cen system; 
e.g. IGuenther fc Demarque! 2000; EgEcnbcrgcr 2004), 10 Gyr for 
Barnard's star, and 0.5 — 3 Gyr for the "thin", 3 — 7 Gyr fo r th e 
"mid", and 7—12 Gyr for the "thick" MEarth samples (see !|4.H . 
These age ranges for the MEarth samples are indicated by the 
lengths of the horizontal lines displaying the percentiles, and the 
MEarth data are shaded gray. The solid curves show rotational 
evolution models for 0.25 Mq stars, fit to the percentiles, Proxima 
and Barnard's star, with the upper curve for the rapid rotators 
(with parameters t^[ sc f ast and f ast ) and the lower curve for 
the slow rotators (parameters T<ji sc s i ow and K w s i ow ). The dashed 
curve shows the result for the rapid rotators if the wind parameter 
(K-w fast ) i s assumed to be the same as for the slow rotators rather 
than allowing it to vary. The dotted curve shows the break-up 
limit from Eq. JBl . 

ing that in reality, there is a large spread in age for each 
of the bins. 

In order to model the evolut i on, w e follow the method 
of iBouvier. Forestini fc Allainl ([1997). We adopt a stan- 
dard wind prescription commonly used for solar-type 
stars (e.g.. IChabover. Demarque fc Pinsonneaultl 119951 : 
Bar nes fc Sofial I1996D. using a modified version of the 
Kawalerl (|1988I) wind law with a = 1 and n = 3/2, intro- 
ducing saturation of the angular velocity dependence at 
a critical angular velocity w sa t: 

!/ \ 1/2 / \ -1/2 

(3) 

By coupling this wind prescription to a stellar evo- 
lution model , here the Lyon NextGen models of 
iBaraffe et all <fl998h . we model the full evolution of the 
angular velocity w as a function of time. This model has 
several free parameters. We assume an initial angular 
velocity u>o, here fixed to the observed velocities in the 
ONC, and this is held constant for a time Tdi sc to ac- 
count for disc-related angular momentum losses ("disc 
locking" ) on the earliest parts of the pre-main-sequence. 
The slow and fast rotator branches are allowed to assume 
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different values of both of these parameters, with the ex- 
pectation that the fastest rotators should decouple from 
their discs earlier than the slower rotators. For the fully- 
convective objects modeled here, we treat the stars as 
solid bodies. The subsequent evolution is then governed 
only by the two parameters K w and w sa t, which control 
the normalization of the wind losses and the saturation 
rate, respectively. 

For a solar-type star, K w can be calibrated by forc- 
ing the model to reproduce the observed rotation rate 
of the Sun. For M-dwarfs there is as yet no analogous 
object with an old, well-determined age and rotation pe- 
riod (although Proxima comes closest to providing it), 
and indeed the rotation rates of the slow and fast rota- 
tor branches appear to not yet have converged even by 
the oldest ages we consider, so instead we simply fit for 
this parameter. 

Saturation is thought to occur when the rotation pe- 
riod beco mes much smaller than the convective overturn 
time (e.g.. lKrishnamurthi et afl H997). These authors in- 
troduced a scaling of w sa t with mass, of the form: 

^sat=W sa t ! Q— (4) 

where r c is the convective overturn timescale at 200 Myr. 
The convective overturn time lengthens for lower mass 
stars, resulting in w sat decreasing with decreasing stellar 
mass (corresponding to slower rotation). 

Typical values of oj sat adopted by other authors for 
the masses we consider here (0.25 Mg) are of order a 
few times the (present-day) solar ang ular velocity (e.g., 
I Sills. Pinsonneault fc Terndrud [2000h . corresponding to 
periods of a few to ten days. 

It is important to note that the convective over- 
turn timescales are no t well-known for fully-convective 
stars; ISills et aLl (120001) extrapolat ed the theoretical val- 
ues from iKim fe Demarquel (|1996f ) to these masses, and 
found that the resulting values of w sa t caused too much 
angular momentum to be lost compared to the obser- 
vations in the Hyades, which were better reproduced 
with a smaller value of u sa t- More r ecent studies of 
empi r ical convective overturn timescales (jPizzolato et al.1 
I2003t iKiraea fc Stepiehl I2007D indicate that r c may in- 
crease sharply around the fully-convective boundary, 
meaning such an extrapolation might under-predict r c , 
and thus overpredict w sa t, as found. It is also possi- 
ble that the Rossby number is no longer the dominant 
factor in determining the operation of the magnetic dy- 
nam o driving stellar winds in fully-convective ob jects 
(e.g. iMohantv fe Basrill2003t [Reiners fe Basrill2007D . 

The uncertainties associated with the Rossby scaling 
have led us to consider several possib ilities for a; s a t. W e 
first consider a value close to that of ISills et al.1 (2000). 
We then allow essentially free adjustment of w sa t f ^5.2.2j) . 
and finally consider a third possibility of there being no 
un-saturated regime, for comparison, in ^5.2.31 

5.2.1. High w sa t model 

If we start by enforcing values of w sa t of a few ujq , the 
model shown in Figure Q3] results. We have performed 
two fits, one enforcing the same value of K w for both 
the fast and slow rotator branches, as would be expected 
in the saturation formalism for solar-type stars, and one 
allowing different values of K w . We consider first the 



former model, enforcing the same value of K w for both 
branches. 

This model is able to reproduce the cluster slow rota- 
tors fairly well, and the oldest ages (Proxima, Barnard's 
star, and the "thick" MEarth sample). The disagree- 
ment with the Praesepe sample and the "thin" sample is 
more concerning, although we note that the age range of 
the latter in practice means it could contain objects as 
young as the Hyades/Praesepe, and that the "thick" and 
"mid" samples could contain objects as young as a few 
Gyr as well. This model does not reproduce the rapid 
spin-down at late times discussed in §5.11 

In the case of the Hyades/Praesepe this model would 
predict that there should be objects rotating as slowly as 
the sun by this age. There are to our knowledge no such 
objects reported in the literature at the present time, but 
we suspect the majority of cluster studies carried out at 
these masses and a ges are not sensiti v e to p eriods this 
long. For example. lScholz fe Eisloffell (|2007H state that 
their sensitivity to periods longer than 1 week is limited. 

The situation for the rapid rotators is much less sat- 
isfactory, with difficulties reproducing the data during 
both the early-PMS and the main sequence phases. We 
first examine the latter. 

As for solar-type stars, the wind prescription we have 
used produces a rapid convergence in rotation rates at 
late times. This does not seem to be the case for M- 
dwarfs. If we instead allow a different value of K w for 
the fast rotators, corresponding to allowing the normal- 
ization of the wind law to be rotation rate dependent, it is 
possible to much better reproduce the observations from 
the age of NGC 2516 onwards. The value of K w found 
here is two orders of magnitude smaller than that for 
the slow rotator branch: these objects experience much 
weaker winds in our model. 

The theoretical interpretation of such a change likely 
lies in the dynamo mechanism operating in these objects, 
as mentioned in ^TJ Of the possibilities mentioned there, 
the a 2 dynamo would possess a strong rotation rate de- 
pendence, for example. 

We now return to the question of the early-PMS evo- 
lution of the rapid rotators. Figure [14] shows that, con- 
trary to the slow rotators, these objects appear to spin 
up more slowly than would predicted simply from con- 
traction. Since the stellar evolution model we are using 
has no rotation rate dependence, our model would only 
be able to reproduce such an evolution by losing angu- 
lar momentum at early-times, in this case by choosing 
a very long disc coupling time of 16 Myr. This seems 
physically unreasonable (and still does not fit the data), 
especially recalling that the rapid rotators are thought 
to be produced from stars which uncoupled from their 
discs early. 

It is possible the prescription for disc locking we have 
used, of a constant angular velocity, is not correct, and a 
more gradual angular momentum removal operates. This 
could result in reduced (rather than eliminated) spin-up. 
The slow rotators do not seem to support this, showing 
apparently relatively little evolution as measured by the 
percentiles out to ages of 5 Myr. It is however impor- 
tant to note that the slow rotation end of the samples 
is subject to contamination from field objects, so it is 
possible this conclusion is a result of the higher level of 
field contamination in the NGC 2362 sample, compared 
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Fig. 15. — As Figure [T4l except for low values of oj sa t- 

to the ONC and NGC 2264. 

Another possible clue may be the proxim i ty of these 
objects to the break-up limit. iHerbst et al.l (|2001l) give 
the following formula for the period corresponding to 
break-up: 



■Pbreak = 0.116 days 



(^/R Q ) 3/2 

(M/M Q ) 1 /2 



(5) 



where M and R are the stellar mass and radius, respec- 
tively. This is plotted as t he do tted line on Figure fL4l 
As stated bv lHerbst et al.l ((2001), many of these objects 
are indeed rotating very close to the break-up rate, with 
many approaching 30% of break-up. This appears to 
persist to ages of approximately 10 Myr. 

In reality, the proximity of these objects to the break- 
up limit means the non-rotating stellar evolution model 
we have used may no longer be applicable, because it 
implies centrifugal forces are likely to play an important 
role in determining the stellar structure. It may there- 
fore be necessary to invoke a more complete treatment 
of stellar evolution, incorporating rotation, to reproduce 
the beh avior of these obj ects, for example the treatmen t 
used by lSills et all (|2000l ) and lDenissenkov et al.l (|2010h . 

5.2.2. Low cj sa t model 



If we relax the assumption made in i)5.2.1l of a value 
of oj sa t of a few w Q , it is possible to obtain an alternate 
solution which better-fits the behavior at late times, at 
the cost of worse agreement with the open cluster data 
around the age of the Pleiades. We show this solution in 
Figure [T5I which has w sa t = 0.65 w©. 

This model appears more satisfactory overall, espe- 
cially considering the age ranges in the field samples, and 
the strong mass-dependence and possible field contami- 
nation in the open cluster samples, meaning some scatter 
for the slow rotators in such a wide mass bin is not en- 
tirely unexpected. The values of the wind parameter K w 
still differ for the fast and slow rotator branches, but now 
by only one order of magnitude. The rapid spin-down at 
late times discussed in £15. II is now produced. 



5.2.3. Zero oj sa t model 

iReiners fc Bas ri (2008) have offered an interpretation 
of the available v sin i data in terms of a saturated wind 
braking law, with a temperature dependence yielding 
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Fig. 16. — As Figure [141 except assuming all objects remain in 
the saturated regime throughout their evolution. 

weaker angular momentum losses for lower-mass objects. 
In their analysis, they assume all the objects rotate faster 
than the saturation rate, corresponding to elimination of 



the dJ/dt 



branch in Eq. ©. This could also 



correspond, for example, to a different magnetic dynamo 
than the solar-type objects have at low rotation rates. 

After the object arrives on the main sequence, the mo- 
ment of inertia and temperature are approximately con- 
stant for a given star, so this angular momentum loss 
prescription can be trivially integrated to show that the 
angular velocity is an exponential function of time. This 
drives the object very quickly to extremely low rotation 
rate. In the context of Figure Q3J where both axes are 
logarithmic, this law corresponds to an exponential de- 
crease in logw versus logt. 

Nonetheless, given the limited quantity of data avail- 
able it is possible to obtain a somewhat satisfactory fit 
using this model, as shown in Figure [TBI 

5.2.4. Overall comments 

The first two models we have considered in this section 
both exhibit a non-saturated spin-down at late-times on 
the slow rotator branch, with the periods evolving follow- 
ing a t 1 / 2 law as for solar-type stars. While we have es- 
sentially enforced the solar-type law by our choice of wind 
prescription, it is important to note that, in reality, dif- 
ferent power laws may be possible depending on the ex- 
act magnetic geometry. Nonetheless this argues in favor 
of a phase of "non-saturated" angular momentum losses. 
Some observational evidence in favor of fully-convective 
M-dwarfs being able to support a non-saturated dynamo 
exists in th e form of a possible rotation-activity corre- 
lation (e.g.. lDelfosse et a l. 1998; Mo hantv fc Basrill2003l: 
although see also lWest fc Basrill2009D . 

We note that t he period of » 1 30 da ys measured for 
Barnard's star bv lBenedict et al.l (|1998|) may also argue 
in favor of a t 1 / 2 or similarly slow spin-down at late times, 
and thus against the model in ^5 .2.31 Al though the ev- 
idence for such a periodicity is stated by iBenedict et al.l 
(1998) themselves as being weak, the MEarth sample 
contains similarly long periods. Barnard's star appears 
to be a thick disc object, and thus i s likely to have an 
age o f approximately 7— 13 Gyr (e.g.. lFeltzing fc Bensbvl 
2008). If we assume Proxima and Barnard's star can 
be treated as evolutionary analogs, the ratio of rota- 
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tion periods of 0.64 is remarkably close to the predicted 
l/y/2 ~ 0.71 from a t 1 ^ 2 law if Barnard's star was twice 
as old as Proxima. The exponential spin-down in Figure 
1161 would produce a much larger change in period. 

Finally, the "pile-up" of objects at long periods seen in 
the MEarth data is difficult to reproduce with an expo- 
nential time-dependence in the late-time evolution. The 
convergence in rotation rates produced by the transition 
to a t 1 / 2 type spin-down law more naturally produces 
this feature, whereas the exponential spin-down instead 
produces a divergence in periods at late-times when we 
use the wind parameter variation with rotation rate in 

Given these lines of reasoning, we believe the model 
in 8 35.2.31 is unlikely, and further has the somewhat un- 
satisfactory feature of producing extreme slow rotators 
( > 1000 days) at th e oldest ages (see also the discussion 
in lReiners fc Basrill20"08l ) . 

Of the remaining two models, we suggest that the ob- 
served spi n-down at late-ti mes (see £|5.1[) and the activity 
lifetimes (West et al. 2008) support the model presented 
in §5.2.21 A comparison of Figures [14] and [15] indicates 
that the largest differences between these two possibili- 
ties (noting that the age range and small number statis- 
tics for the "thin" sample does not strongly constrain the 
model at present) are manifest at the age of the Hyades 
and/or Praesepe, with one model predicting the exis- 
tence of objects rotating as slowly as the sun, and the 
other being consistent with the presently-available ob- 
servations. More generally, the slowest rotation period 
seen at this age constrains the value of cj sat . We believe 
these arguments motivate extending the observations to 
search for slowly rotating M-dwarfs in these clusters. 

6. CONCLUSIONS 

We have presented a sample of rotation period mea- 
surements for 41 field M-dwarfs inferred to lie below the 
fully-convective boundary, from the MEarth transit sur- 
vey. We measure a wide range of periods, from 0.28 days 
to 154 days, where the latter is commensurate with the 
sensitivity limit resulting from the typical baseline of 
the observations, finding that the kinematically young 
(thin disk) objects rotate faster than the kinematically 
old (thick disk) objects. 

Combining the available measurements from the liter- 
ature with the present sample shows two interesting fea- 
tures in the mass dependence of rotation for M-dwarfs at 
the oldest field ages. One is the familiar discontinuity at 
approximately 0.6 Mq seen in many open cluster stud- 
ies, where the mass-dependence of the rotation periods 
steepens below this mass. Contrary to the findings of the 
open cluster studies, the field stars below 0.6 M Q appear 
to rotate slower with decreasing mass. 

The second feature is the onset (moving down in mass) 
of rapid rotation at old ages for some objects below 
0.35 Mq, which is also seen in the vsmi surveys. Our 
rotation periods confirm that not all such field M-dwarfs 
are rapid rotators, and there exists a significant, or even 
dominant, population of slow rotators that appears be- 
tween the ages corresponding to our thin and thick disk 
samples. In addition to this confirmation of the results 
of vsini studies (e.g. iDelfosse et al.lll"998l ). we can now 
estimate the rotation rates for the slow rotators, finding 
them to rotate approximately two orders of magnitude 



slower than the remaining rapidly rotating objects at the 
age corresponding to our thick disk sample. There is no 
discontinuity in the rotation rates crossing 0.35 M for 
the slow rotator branch at the oldest ages. 

Comparing the results for the partially-radiative and 
fully-convective objects indicates that the spin-down 
time increas es moving to lower masses . The "young disk" 
sample from iKiraga fc Stepienl (|2007l ) has a comparable 
fraction of rapid rotators to the "mid" and "thick" sam- 
ples from the present survey, indicating that the spin- 
down times for these populations are comparable to the 
respective sample ages: the "young disk" sample age for 
the partly-radiative stars (a few Gyr) and the thick / 
"old disk" (probably > 5 Gyr) for the fully-convective 
stars. 

We have explored wind models and were able to re- 
produce these features by adopti ng a very low satura- 
tion velocity within a standard (|Chabover et al.l 119951 : 
iBarnes fc Sofial[r996t l wind prescription, although it was 
necessary to assume a rather large rotation rate depen- 
dence of the wind parameter. In reality, this additional 
rotation rate dependence probably corresponds to mod- 
ifying the power law exponent n in the wind model, and 
thus the magnetic field geometry, although it is beyond 
the scope of this paper to speculate on such a change. 
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0.0082 


Thick 


LHS2088 


08 h 59 m 56 s 


1 


+72°57'36" 


1999.3 


0.973 


-0.035 


0.0726 


0.0034 


1 


8.99 


0.15 


0.18 


138 


0.0033 


Thick 


LHS267 


09 h 20 m 57 s 


9 


+03°22'06" 


2000.1 


0.316 


-1.139 


0.0608 


0.0041 


1 


8.54 


0.23 


0.25 


118 


0.0027 


Thick 


LHS2175 


09' l 42 m 23 s 


3 


+55°59'02" 


1999.0 


-0.703 


-0.517 


0.0651 


0.0087 


1 


7.55 


0.31 


0.33 


72.8 


0.0069 


Thick 


GJ403 


10 h 52 m 04 s 


4 


+ 13°59'51" 


1998.0 


-1.115 


0.198 


0.0839 


0.0257 


1 


7.81 


0.25 


0.27 


96.8 


0.0038 


Thick 


G119-62 


ll h ll m 51 s 


8 


+33°32'11" 


1998.4 


-0.164 


0.110 


0.0683 


0.0106 


1 


7.51 


0.31 


0.33 


7.79 


0.0197 


Mid 


LHS316 


ll h 50 m 57 a 


9 


+48°22'40" 


1999.0 


-1.540 


-0.960 


0.1221 


0.0029 


3 


7.66 


0.17 


0.19 


132 


0.0042 


Thick 


LHS6220 


12 h 05 m 29 s 


7 


+69°32'23" 


1999.2 


-0.447 


-0.040 


0.0602 


0.0134 


1 


7.91 


0.24 


0.26 


96.2 


0.0064 


Thick 


GJ1156 


12' i 18 m 59" 


4 


+ 11°07'34" 


2000.2 


-1.263 


0.199 


0.1529 


0.0030 


3 


7.59 


0.14 


0.16 


0.491 


0.0142 


Thick 


GJ493.1 


13 h 00 m 33 s 


5 


+05°41'08" 


2000.2 


-0.934 


0.215 


0.1231 


0.0035 


3 


7.68 


0.17 


0.19 


0.600 


0.0028 


Thick 


LHS2995 


14 h 53 m 37 s 


2 


+ 11°34'12" 


2000.3 


0.076 


-0.738 


0.0535 


0.0041 


1 


9.67 


0.15 


0.17 


88.3 


0.0074 


Thick 


GJ1204 


16 h 36 m 05 s 


6 


+08°48'49" 


2000.4 


-0.519 


-0.151 


0.0652 


0.0042 


1 


8.53 


0.21 


0.23 


6.41 


0.0033 


Thick 


NLTT43307 


16 h 40 m 06 s 





+00°42'19" 


2000.3 


0.174 


-0.163 


0.0890 


0.0023 


1 


8.23 


0.18 


0.20 


0.311 


0.0043 


Thin 


LHS6311 


16' l 40 m 20 s 


7 


+67°36'05" 


1999.4 


-0.262 


0.361 


0.0747 


0.0043 


1 


8.97 


0.15 


0.17 


0.378 


0.0091 


Thick 


G141-29 


18 h 42 m 45 s 





+ 13°54'17" 


1999.4 


-0.025 


0.347 


0.0933 


0.0115 


1 


7.57 


0.15 


0.17 


8.07 


0.0132 


Thick 


LHS3461 


19' l 24 m 16 s 


3 


+75°33'12" 


2000.4 


0.374 


0.591 


0.0903 


0.0051 


1 


8.98 


0.12 


0.15 


114 


0.0121 


Thick 


GJ1243 


19 h 51 m 09 s 


3 


+46°29'00" 


1998.4 


0.188 


0.266 


0.0841 


0.0024 


1 


7.79 


0.24 


0.26 


0.593 


0.0085 


Mid 


GJ791.2AB" 


20 ,l 29 m 48 s 


4 


+09°41'20" 


2000.4 


0.665 


0.132 


0.1138 


0.0019 


3 


7.33 


0.21 


0.23 


0.346 


0.0038 


Mid 


GJ856A 


22 h 23 m 29 s 





+32°27'33" 


1998.5 


0.255 


-0.208 


0.0622 


0.0100 


2 


6.07 


0.25 


0.27 


0.854 


0.0070 


Thin 


LHS543a 


23 ,l 25 m 40 s 


2 


+53°08'06" 


1998.9 


0.986 


0.328 


0.0404 


0.0031 


1 


9.00 


0.28 


0.30 


24.8 


0.0077 


Thick 


LHS3978 


23' l 35 m 41 s 


3 


+06°11'21" 


2000.7 


0.566 


0.268 


0.0416 


0.0032 


1 


10.19 


0.15 


0.18 


1.67 


0.0055 


Thick 



TABLE 1 

Rotation periods, kinematics and stellar parameters for the MEarth 

SAMPLE. 



a Coor dinates from 2MASS, ep och of measurement as specified in the "epoch" column. 
b From lLepine fc Sharal ([200l . 

c Reference for parallax measurement: (1) Yale Parall ax Catalog, van Altena et aTl (Il995[); (2) The H ipparcos catalog. IPerrvman et al.l (|1997l) ; (3) NASA 
NStars database ( |http://nstars.arc.nasa.go"v7|; (4) iDahn et al.l (j2002t ) . Compiled bv lLepind (|2005l ). 

d From 2MASS, converted to the CIT system of lElias et al.1 (|1982L Il983[) using th e transformation from the 2MASS explanatory supplement. 

c Using the -ftf-band absolute magnitude and the relation of iDelfosse et al.l (|2000D . 

f Derived from M using a polynomial fit to the empirical mass-radius data of[Ribas (2006)). 

s While we have not attempted to estimate formal period uncertainties, our simulations indicate that the uncertainty in the periods is < 1% for P < 10 d, 
« 1% for 10 < P < 20 days, 2% for 20 < P < 50 days, 5 - 10% for 50 < P < 100 days, and 20 - 30% for P > 100 days. 
h Semi-amplitude of modulation seen in the MEarth bandpass. 
1 Kinematic population, see ^4.1 1 

J Visual binary, unresolved in these observations; see iBenedict et "all (|2000D . 



